Nonperturbative master equation solution of central spin dephasing dynamics 
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We solve the long-standing central spin problem for a general set of inhomogeneous bath couplings 
and a large class of initial bath states. We compute the time evolution of the coherence of a central 
spin coupled to a spin bath by resumming all orders of the time-convolutionless master equation, 
thus avoiding the need to assume weak coupling to the bath. The fully quantum, non-Markovian 
solution is obtained in the large-bath limit and is valid up to a timescale set by the largest coupling 
constant. Our result captures the full decoherence of an electron spin qubit coupled to a nuclear spin 
bath in a GaAs quantum dot for experimentally relevant parameters. In addition, our solution is 
quite compact and can readily be used to make quantitative predictions for the decoherence process 
and to guide the design of nuclear state preparation protocols. 



Since the central spin model was first introduced by 
Gaudin [l[ several decades ago, it has appeared in diverse 
physical settings such as spin dynamics in disordered in- 
sulators [2|, interacting electrons in metallic grains |3|, 
and more recently in semiconductor spin qubits |4l4l3|. 
It has also been widely studied in the context of inte- 
grable systems, where close connections to BCS theory 
and related theories of pairing forces have been discov- 
ered [3|, [ij]. In addition, several variants of the central 
spin model have served as toy models for compa ring and 



in 



contrasting various master equation formalisms [15- 

This broad interest in the central spin model has led 
to a plethora of disparate approaches to computing its 
dynamics. Since its inception, it has been recognized as 
an integrable system which admits a solution via Bcthc 
ansatz [l|, [j, Q- However, this yields only a very compli- 
cated ground state wavefunction, and progress in extract- 
ing dynamical information for the central spin has been 
limited; to date only modest bath sizes of up to 30 bath 
spins have been treated in this manner [8| . Sophisticated 
numerical recipes for calculatin g th e central spin evolu- 
tion have been developed [Sl, lisl, |l9| , but success tends to 
require either small baths (tens of spins) or the assump- 
tion that the initial bath state is completely unpolarizcd. 
These requirements are too restrictive for many applica- 
tions in the context of III-V semiconductor spin qubits. 



where the number of nuclear spins ranges from 10 to 
10^, and polarized baths are employed to facilitate qubit 
operations and extend coherence times 20|, [21 1 . 



This state of affairs led to the development of several 
approximate analytical methods to compute the central 
spin evolution. Some of these methods employ an effec- 
tive pure-dephasing Hamiltonian which is derived per- 
turbatively from the central spin Hamiltonian using a 
canonical transformation [7|, [ol llO| . This approach cul- 
minated with a nonperturbative solution of the effective 
Hamiltonian dynamics [9|, llOf and has been successful in 
describing spin echo experiments in the case of unpolar- 
ized baths [22|. However, an uncontrolled approximation 



in the effective-Hamiltonian derivation makes it unclear 
when this approach is valid, and it has yet to be ex- 
tended to more general bath states. Methods employing 
generalized master equations, on the other hand, offer a 
controlled approximation and naturally describe polar- 
ized bath states, but only perturbative treatments have 
been given so far [y, |ll|, |l3|, [IGI, [l? | , leading to solutions 
which arc valid only outside the regime relevant for many 
semiconductor spin qubit experiments J22h25i |. 

In this Letter, wc solve the central spin problem us- 
ing the time-convolutionless (TCL) master equation for 
a general set of inhomogeneous coupling constants and a 
large set of initial bath states, including both polarized 
and unpolarizcd baths. The TCL equation is an exact 
equation for the reduced density matrix of a system cou- 
pled to a bath; although this equation is time-local, it 
incorporates the full bath dynamics [ly, ujl, l26[ . With 
only a very modest condition on the magnetic field, we 
give a closed-form solution describing the evolution up 
to a timescale set by the largest bath coupling. We are 
therefore presenting an exact solution to the central spin 
model as it pertains to gated GaAs spin qubits since this 
temporal window contains the entire decay of the electron 
spin coherence in the low magnetic field regime where 
the central spin model gives a good description of the 
physics 27|. Although we will focus on the example of 
spin qubits, our results are potentially applicable to any 
Gaudin-type central spin problem. This result is also im- 
portant to the general study of open quantum systems; 
wc are not aware of other examples involving a large, 
nontrivial and highly non-Markovian bath where an all- 
orders rcsummation of a master equation expansion is 
performed. 

The central spin model is comprised of a central spin 
coupled to a spin bath via a Hciscnberg interaction. As- 
suming a nonzero external magnetic field, the Hamilto- 
nian is 



H = J2 ^kS ■ h + ns' + J2 "^kik 



(1) 



with S denoting the central spin operator, I^ the bath 
spins, Ak the (hyperfine) couphngs, and fl and ujk the 
central spin and bath spin Zeeman energies, respectively. 
In the context of semiconductor electron spin qubits, the 
Ak are determined by the shape of the electron wave- 
function envelope, but we will leave the Ak completely 
arbitrary. We refer to J^k ^kS^I^ as the Ovcrhauser 
term, and Vs=^ Y.k ^k{S^Ik + •5'"-^^!") as the flip-flop 
term, where S^ = S*^ ± iS^ , and similarly for I^ . The 
total interaction energy is A = ^^. Ak, and the number of 
bath spins appreciably interacting with the central spin is 
A^ = -4^/(X)fe ^D- Roughly speaking, the bath produces 
an effective (Ovcrhauser) magnetic field, the magnitude 
of which is controlled by A and the bath polarization, 
and about which the central spin precesses, while the Ak 
set the scale for the precession of individual bath spins 
about the central spin. 

The TCL equation is an exact equation for the time 
evolution of the reduced density matrix of a system cou- 
pled to a bath [26[. Although this equation contains full 
memory of the bath dynamics, unlike equations such as 
the Nakajima-Zwanzig equation [a, |ll|, |l3[ , it has the at- 
tractive feature that it is a time-local ordinary differential 
equation. Working in an interaction picture defined with 
respect to Hq = J^k^kS^Ik + ^S'^ + J^k^klk and de- 
noting the total density matrix in the interaction picture 
by p(i), the TCL equation has the form 



-pp{t) = j2ic,,{t)pp{t). 



dt 



(2) 



The operator P projects the full density matrix onto the 
reduced density matrix of the system; its precise defini- 
tion will be given shortly. Eq. ([2]) is defined in terms of 
a perturbative expansion in T4f which is nominally con- 
trolled by the quantity A/^. Ultimately our solution 
will not require this expansion to be convergent since we 
will sum the entire series, so A/fl need not be small. The 
nth-order kernel A^„ (t) encapsulates full bath effects aris- 
ing from nth-order fiip-flop processes in which the central 
spin fiips n times with one or more bath spins. It can be 
expressed as an integral of ordered cumulants involving P 
and the interaction-picture Liouville operator L (defined 
by p = —iLp); for example the second-order kernel is 
IC2{t) = -J*dt'PL{t)L{t')P. The rules for constructing 



26| 



the higher order kernels can be found in Ref. 

We will assume that the initial density matrix sepa- 
rates into system and bath components, p{0) — ps{0) <^ 
Pr(0), and we will use a set of "correlated projectors" 
[l6| , in which case the action of P on a matrix M is PAI ~ 
J2a TrslnoA/} ® j^IIq, where the Ua are a set of bath 
projectors satisfying IlaIlp = Sap^p and J2a^oi = ^, and 
A/'a = TrsjIItj}. The 11^ allow us to write the reduced 
density matrix of the system as a sum of independent 
degrees of freedom: ps = J2a TrsjIIap} = J2a Ps- The 



which in turn implies that the choice of the IIq, will con- 
strain the possible initial bath states. An appropriate 
choice of the 11,-, can either simplify or vastly improve the 
convergence of the TCL equation depending on the sym- 
metries of the Hamiltonian and initial bath state J13l Il6l | ; 
for now we leave the IIq, completely general. 

Since we want to compute the off-diagonal component 
of ps{t) (coherence function), we multiply Eq. ^ by 5+ 
and trace over both system and bath to obtain 

oo 

ps,-+it) = Ti-s{S+psit)} = ^Tr{5+/C„(t)Pp(t)}. 

(3) 
The summand on the right-hand side of this equation is 
comprised of integrals of terms with the general struc- 
ture Ti{S~^CiPC2P ■ ■ ■ ^rPp{t)}, where £i represents a 
string of Liouville operators L(ti-^)L{ti^) .... If one as- 
sumes that Ua is such that bath correlators of the type 
Tr^jIlQ,/^ I^ . . .} vanish unless they contain equal num- 
bers of raising and lowering operators, then it is straight- 
forward to show [29[ that in the case of the central spin 
model H]), such terms factorize: 

Tt{S+CiPC2P . . . CrPpit)} 

= E Tf^Pt-+it) n Tr{5+A5-n„}. (4) 



i=l 



Since every term on the right-hand side of Eq. ([TT|) fac- 
torizes in this way, we can expand the left-hand side as 

Ps, h(0 = X^a Ps-+i^) ^^'^ separately equate each term 

of the a-sum. The resulting set of equations is readily 
solved: 



p£+(o=P5':i+(o)cxpK]e("no 



(5) 



where gi''\t) = /^ dt'Tr{5+/C„(t')^-n„}. 

To calculate Gn (t) we must first compute correlators 
of the type Tr{S'+£S'~nQ,} where £ is an arbitrarily long 
string of Liouville operators. Restric ting ourselves for 
simplicity to the case w^ = w, we find [29| 

Tr{5+L(4j . ..LihJS-Ua} « ^e-^''- I^l.i-^)''K 



4? 



fc=0 



J2[i)T^^B{ih+h-fiiaih-h+y-'^}. 



(6) 



We have defined the operators h^ = J2e -^i-^l and fre- 
quency J7q, = il — w + Tv B{h''Ua\ / J^cn and we have as- 
sumed that [ft,^,!!^] = 0. The latter condition ensures 
that the correlator vanishes for odd numbers of Liouville 
operators, which in turn implies that Qn vanishes for 
odd n. There are two approximations being made in 
Eq. (HH. The first approximation assumes t <C 1/Amax, 
form of the TCL equation we use requires Pp(0) = p(0), where Amax^A/N is the largest coupling, and this leads 



to the timc-depcndcnce appearing only as a phase factor 
in Eq. ^Q. For spin qubits in GaAs with TV = lO*', 
l/Ajnax can be on the order of 10/iS, long enough to 
capture the full decay of the electron spin coherence for 
exp erimentally relevant values of the magnetic field [22|- 
25| . For comparison, we can also consider spin qubits in 
Si [30[, in which case l/A^ax is on the order of 250/iS. 

If this were the only approximation, then the right- 
hand side of Eq. (f2T|) would have an additional sum 
over permutations of the tb^ , but only a certain subset 
of these permutations were kept in Eq. (pTj) . Retaining 
only this subset amounts to keeping the leading order 
terms in the flat ^ 1 limit at each order of the TCL 
expansion. To illustrate the nature of this RPA-like ap- 
proximation, we consider its effect on the lowest-order 
terms. If we kept all permutations, then the second-order 
term would have the form G2 v-) ~ ^a'^i^j ^a) while the 



cases in (|2T 
where Sa = 



would have the form Q^"'{t) ^ i5^D(i, fie) 



and u(i,fia) = i^at + 6"**^°* - 1, 
iiQat -{6 + 2zr2„t)e-*f^°* + 6. The 



RPA-like approximation amounts to taking u — > iflat 
and V — >■ — fJ^t^. Self-consistency of the approximation 
requires (5^ <C 1, which imposes a lower bound on the 
magnetic field (corresponding to a few mT for gated dots 
in GaAs). The necessity for 6a <S^ 1 can be seen from u 
and V by noting that we can only neglect the linear (in 
t) term in v if it is small compared to the linear term 
in u. Later on, we will see that this approximation cap- 
tures the envelope of the coherence function, and we will 
also find that we can relax this approximation to a large 
degree by keeping the full form of Q2 (t)- 

With Eq. ([21]) in hand, it is straightforward to assem- 
ble these correlators into the function Qn (t) using the 



fourth-order terms stemming from the q — 2, k ~ 0,2 rules for constructing the TCL ordered cumulants [2£ 
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i=l fe=0 



HE rnTrsiih^h-ruaih-h+r-'^}. 



(7) 



In this expression, {qi} € V{n/2) means that {qi} is an 
ordered integer partition of n/2, with r being the num- 
ber of Qi comprising the partition [3l| . To evaluate this 
expression, we make perhaps the simplest choice for the 
projectors: Ua=Il^ = \x) {x\, where |x) = <S>k \^k,rnl) 



is a product of eigenstates of the I^ {Ik {Ik 



1) and m^ 



are eigenvalues of I^ and /|). In this case, A/'q = 1. This 
particular choice is well suited to applications pertaining 
to spin qubits [6|,|ll|,ll3|. With an initial bath density ma- 
trix of the form p_b(0) = X^v Pxx \x) (xl; it is possible to 
perform the various sums in iJ23|) despite their complexity 
[29J , and we find to leading order in the limit of large N 

that the coherence function W{t) = ps, ^{t)/ps, i-(0) in 

the Schrodinger picture is 



W{t) 



E 

X 



p^,(d+-ci-)e^("+''x)* 



d+e 



m;idt-d^) 



d^"''^ 



n7(4-dx)' 



(8) 



where h^ = Trs{n^/i^} is the Overhauser field as- 
sociated with the state |x), d^ = (xl ^^^^ Ix) — 
J2i A'j [Ig{I(i + 1) — mf{m^ ± 1)] quantifies transverse fluc- 
tuations of this field, and Vt^ = fi— w-t-Zif^ is the difference 
between the effective Zeeman energy of the central spin 
{V, + h^) and the bath spin Zeeman energy (w). Equa- 
tion (|40p is the main result of this paper; it describes the 
envelope of the coherence function of the central spin for 
an arbitrary set of couplings Ak and for a large set of 
initial bath states up to time t < ^/Amax- 

As a first example, we consider the case of a uniformly 



polarized bath (32| with all bath spins having the same 
total angular momentum /, and for which /if, is the same 
for all X liEI- Writing /ij = Alp and np^il-u + AIp, 
where p = -^ ^ p^^ J^k "^fe i^ ^^^ average polarization 
of the bath, with p = denoting an unpolarized bath 
and |p| = 1 maximal polarization, we find 

„i{n+Aip)t 
WW = 77^^ Z 777-T^ (9) 



pcos 



f2Ipt\ 



ipli 



m 



with Pi ^ / + 1 - ir Ex Pxx Y.k{^4? (for / = 1/2, 

p\ = I), and Tp = ANVtp/A^. Setting p=0 in (0) yields 
the zero-polarization result, W{€) = 1/(1 — 2iIp]j./Tp), 
obtained in previous works [9|, [lOl ll3| using less rigor- 
ous methods. The left panel of Fig. [T] shows that the 
decoherence rate increases with increasing / due to a 
corresponding increase in the number of bath degrees 
of freedom. For electron spin qubits in a GaAs nuclear 
spin bath (/ = 3/2) where the magnetic field is typi- 
cally on the order of 100 mT [^ so that A/Vt « 30 
and To « (2/15)A^/^, it is clear that the coherence de- 
cays almost completely before time t = 1/Amax ~ N/ A 
is reached. The right panel of Fig. [T] depicts the extent 
to which the decoherence time increases with increas- 
ing bath polarization, a well known effect which can be 
understood in terms of a reduction of phase space for 
flip-flops. It is also evident that positive net polarization 
leads to longer decoherence times compared with nega- 
tive polarization, since for the latter the decrease of \^p\ 





FIG. 1: (Color online) Coherence function (Eq. ((9|). Left 



panel: changing bath total spin / with p = and p^ 



(1 + 



1). From top to bottom, / takes the values i, |, |, |, |. Right 
panel: changing bath polarization p with I ~ ^, A/fl — 
30, and pi = |. From top to bottom, p takes the values 

3 j_ 1 _1 J_ ^J_ Q 

4' 2' 4' 4' 10 ' 10 ' 



facilitates virtual flip-flops. Our results provide quantita- 
tive predictions for the enhancement of decoherence time 
resulting from bath polarization. 

A salient feature of the uniform-polarization result, Eq. 
([9]) ; is that it depends on the couplings only through the 
quantities A and N . This implies that any set of cou- 
plings which yield the same values of A and N will give 
rise to the same central spin evolution for t <C N/A. For 
example, this evolution should be reproduced by a model 
in which all the couplings are equal, Ak = A/N, the so- 
called "box" model, which is exactly solvable (see e.g. 
[13| for the solution in the case of a polarized bath). A 
comparison of Eq. (j9|) with the exact box model solution 
is shown in the left panel of Fig. [5J and it is evident that 
the two solutions agree very well. This insensitivity to 
the particular values of the Ak on timescalcs t ^ N/A 
was anticipated in [l3| based on energy-time uncertainty; 
here, we have given a direct proof of this result, and we 
show below that its validity requires uniform polariza- 
tion. The exact box model solution is not known to have 
a closed form, so the fact that Eq. ([9|) constitutes a very 
good approximate closed-form solution is an added bonus 
of the present work (Eq. (|40|) ). 




20 40 60 80 100 120 140 
t/To 



0.2 

t/To 



FIG. 2: (Color online) Left panel: Coherence function from 
Eq. ((SJ (blue, dashed) vs. exact box model solution from [13] 
(red, solid) with / = f, A/Q. = 30, iV = 10*, P= jg- Right 
panel: Zoom-in of left panel with modified coherence function 
W{t) (black, dotted) included as well. 

The right panel of Fig. [5] reveals that Eq. ^ does not 
capture a high-frequency, small-amplitude modulation 



FIG. 3: (Color online) Coherence function with uniform po- 
larization Eq. ((5)) (red, solid) vs. non-uniform polarization 
Eq. gni) (blue, dashed) with / = |, A/Q = 30, and (left 
panel) P = iij, Tp = |^, (right panel) p = |, rp = i|^. 



exhibited by the exact box model solution [13|. However, 
it turns out that it is easy to correct for this by relaxing 
the RPA-like approximation at second order of the TCL 



flip-flop expansion, in which case one finds that Q2 



(x) 



1 
4n 



[if2^,f -|- e 



-iUyt 



l](dj + d^), leading to the modi- 
fied coherence function W = 14^exp[(e'*^p* - l)/{flpTp)]. 
This modified function reproduces very well the small 
modulations as is evident in the right panel of Fig. [2] 



The condition S-y 



<S^ 1 ensures that higher-order 



corrections to these modulations are negligible. 

We stress that the coherence function given in Eq. (|40|) 
can describe more general bath polarization states. To 
illustrate this, we consider a simple example of a non- 
uniformly polarized initial bath state where the bath 
spins possess an average polarization which depends lin- 
early on their coupling: (rnfc) = X^v Pxx"^fc = NIpAk/A. 
This describes a bath configuration where the spins clos- 
est to the central spin are more polarized, while those 
further away are less polarized. In the case of spin 
qubits in quantum dots, this is qualitatively a physi- 
cally plausible configuration since the nuclear spins are 
typically polarized throu gh manipulation of the central 
electron spin [20|, |2l|, |33I436| . This time, the coherence 
function computed from Eq. (j40p depends not only on 
J2k^k and X^fc^fci but also on J2k^h meaning that 
the result now depends on details of the distribution 
of the Ak- Assuming the Ak are distributed in accor- 
dance with a two-dimensional Gaussian wavefunction of 
the central spin electron, and taking the continuum limit 
(e.g. Ek fiAk) ^ /o"-"^ p{A)f{A)dA with p{A) ^ ^ 
and Amax = ^), we obtain 



Wit) = 



pe 



i{n+AIp)t 



pcos 



f8Ipt\ 
\3r, J 



^p\ sin 



(8Ipt\ 
\ 3rp J 



(10) 



Fig. [3] shows that while for lower bath polarizations 
this result is similar to what we found in the uniform- 
polarization case, for larger polarizations, the differences 
between the two solutions become quite pronounced. 
This illustrates how our results can be used to distin- 
guish between different narrowed polarization distribu- 



tions produced using empirical nuclear state preparation 
schemes [20|, |2l|, |3J, |35[ which are not yet understood 
microscopically. 

In conclusion, we have presented a non-pcrturbative 
solution to the central spin problem. Our solution is an 
elegant, closed- form expression which is valid for an ar- 
bitrary set of couplings to the spin bath and for a large 
class of initial bath states. In the context of GaAs spin 
qubits, our result solves completely the pure-hyperfine 
decoherence problem and can be used to quantitatively 
predict decoherence times and to better understand nu- 
clear state preparation protocols. 
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SUPPLEMENTARY INFORMATION 

This supplement shows exphcitly the intermediate algebraic steps which arise at different stages of the derivation 
of the central spin coherence function. 

FACTORIZATION OF BATH CORRELATORS 

In the main text, it is shown that the derivative of the coherence function in the interaction picture can be expressed 
as 

oo 

ps.-+{t) - TTs{S+ps{t)} = J2 Tr{5+/C„(<)Pp(t)}, (11) 

where /C„ (f ) is the nth order TCL kernel. This quantity is given in terms of ordered cumulants involving the Liouville 
operator L(t) and the projection operator P. For example, at fourth order we have 

ICi{t) = f dh r dt2 f ' dt3\PL{t)L{h)L{h)L{t3)P - PL{t)L{h)PL{t2)L{t3)P 



Jo Jo Jo '- 

-PL{t)L{h)PL{h)L{t3)P - PL{t)L{t3)PL{h)L{t2)p] . (12) 

The rules for constructing these kernels can be found in Ref. [l| . For the central spin Hamiltonian, the explicit action 
of the Liouville operator on a matrix p is 

^ fp++ P+-\ ^ 1 fh+-p-+ - P+-h-+ h+^p^^ - p++h+\ .^^ 

\P-+ P--J 2\h^+p++- p—h^+ h^+p+^- p^+h+^J ' 

with 

f. 

The summand on the right-hand side of Eq. (|lip is comprised of integrals of terms with the general structure 

Tv{S^CiPC2P ■ ■ .CrPp{t)}, where £,; represents a string of Liouville operators L{ti^)L{ti^) We will now show 

how such terms factorize by focusing on the case r = 2. Because of the particular form of the central spin Hamiltonian, 

when we expand Ty {S'^ LiP L2P p{t)} using Eqs. ([T^]) and P^ . each resulting term depends on either {PL2PP) y or 

{PC2Pp)^ Terms depending on {PC2PP) h always contain equal numbers of bath raising and lowering operators, 

while terms depending on {PC2Pp)^ contain unequal numbers. If we assume IIq is such that TrsjIIa/^ //^ . . .} 

vanishes for unequal numbers of IJ' and I^ , then only the terms with {PC2PP) t- are non- vanishing, and we have 

TT{S+CiPC2Pp{t)} = TT{S+Ci{PC2Pp)-+S-} = TT{S+cJY.^^B{TlaC2Pp}®^TlcA S-} 

= Y. ;^Tr{5+£i5-n„}Tr{5+/:2Ppn„} = ^ ^pgl+Tr{5'+/:i5-n„}Tr{5+£25"n„}. 

a a ^ 

(15) 
This factorization generalizes straightforwardly to correlators with any number of P insertions: 

Ti-{S+C^PC2P . . . CrPpit)} = J2 TrTP^sl+ii) \[ Tr{S+ C,S-n^}. (16) 

The assumption that TrslIlQ./^ I^ . . .} vanishes for unequal numbers of raising and lowering operators is an assump- 
tion that restricts the possible choices we make for 11^, and thus the possible choices of initial bath state. More 
specifically, if we represent IIq as a matrix in the basis of eigenstates of /i^ , then this restriction forces many of the 
off-diagonal elements to vanish. This restriction is rather mild and still leaves us with a large class of initial bath 
states including those relevant in the case of electron spin decoherence in the presence of a polarized nuclear spin 
bath. 



EVALUATION OF BATH CORRELATORS 



The factorization shown in Eq. (jl6[) ahows us to reduce the problem of computing correlators of the form 
Tr{S^CiPC2P ■ ■ ■ CrPp{t)} to computing correlators of the form TrlS'^CS^Ha}- The latter can be evaluated 
straightforwardly: 

-^Tr{^+L(i')i(ti)...L(t2,-i)5-n„} = 1^^5(")({t,J), (17) 

with 

di^^ita,}) = ^TrB{/l-+(tai)/H-(^a2)---/l--+(ia2fc-i)/H-(W)na/H-(W + l)'l-+(^a2fc+2)---/H-(^a2,-i)/l-+(ia2j} 

(18) 

Here, the sum X]|a > means that we sum over all permutations of the a^, where i = 1 . . .2q and a^ G [0,(7 — 1], such 
that ai < 02 < . . . < a2k and 02^+1 > a2fc+2 > ■ • ■ > 02^. We define to = t'. At this point, we further restrict IIq to 
be such that [IIc/i^] = 0, and we focus on the short-time evolution, t < 1/Amax, in which case we may neglect the 
Ak appearing in the exponents in the expressions for h^ and h |- in Eq. (|14p . In this case, gj^ simplifies to 

g[")({t„J) ^ ^e"'"°(*"i-*''^+--*''2.-*»2.+i+*»2.+2 — ■+*<'2,)Trs {{h+h-)''Ila.{h-h+y-''} , (19) 

where /i^ = J2e -^i^f ■ 

The second (RPA-like) approximation we will make is to keep only those permutations of the a^ which yield the 
dominant terms at order n after the integrations over the <a. are performed. Retaining only this subset of permutations 
amounts to keeping the leading order terms in the Q,at 3> 1 limit at each order of the TCL expansion. For example, 
one such dominant permutation is given by iai = t' , ta- = U-i, ta ~ t2q+2k-ji i = 2 . . . 2k, j = 2k + 1 . . . 2q, which 
leads to the time-dependent factor 



-»f2«E-Ii(-i)"'to 



(20) 



In fact, the nested time integrals (see for instance Eq. ([T2|) ) are such that the time-dependent factors in the dominant 
(in the Jlat 3> 1 limit) terms of T^{S^ CS~Iia] always have the form shown in Eq. ([20|) . This means that we need to 
sum over only those permutations of the ai which preserve Eq. \2Q\ . There are (2^.) different choices of the a^, and 
(^) of these will lead to the same time dependence as in Eq. ([20| . This counting can be understood by noticing that 
the quantity t — ti-\-t2 — ■ ■ ■ — t2k-i — t2k + ■ ■ ■ + t2q is invariant under pairwisc swaps of t2i — t2i+i with t2j ~ t2j+i- 
All such swaps which preserve the condition that the first 2k times in this quantity arc monotonically increasing from 
left to right and the last 2q — 2k times arc monotonically decreasing from left to right should be counted as separate 
contributions to the time-dependent factor. Restricting ourselves for simplicity to the case LOk = ^, wc therefore have 

Tr{5+L(taJL(taJ...L(ta,J5-n4 « Ig-^c E^li (-!)"•*=. ^ (^\TvB{{h+h-fn^{h~h+y-''}. (21) 

As explained in the main text, this RPA-like approximation is controlled by the small quantity 5^ = 'jWcT^ which 
places a lower limit on the external magnetic field. 

Now that wc have obtained Eq. (pij) . wc need to combine these correlators according to Eq. (IT51) and integrate the 
result over the ti. Since the time-dependent factors combine the same way for every term at order n, these integrations 
are easily performed all at once: 

, dt' dh dt2... df„_ie-'^^°(*'-*^+-+*"— *"-) = V4^ 7f • (22) 
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Defining gr'{t) = /^ di'Tr{5+/C„(i')^"n„}, we then find 

JLY^" V 1 (-1)^+1 -^/g. 

4^ J ....^,w,.nU'?.! r^fr 11A.I, 



^i"^w= i^ E Tf^^^^nEr; Tr.{(.-/.-)^n.(.-.-)-n, (23) 



where {qi} € P{n/2) means that {qi} is an ordered integer partition of n/2, with r being the number oi qi comprising 
the partition. For example, ii n = 6, then the possible choices of the g, are {qi} = {3}, {2, 1}, {1, 2}, {1, 1, 1}, 
and r takes the values 1,2,2,3 respectively. Summing over these partitions and including the combinatoric factor 
(n/2)!/(r ]^[^j^ g^!) is tantamount to summing over the different terms appearing in the nth order TCL kernel (as 
dictated by the rules for constructing ordered cumulants [l[), for example the four terms in Eq. P^ . An additional 
factor of (— l)"/2+''+i also comes directly from the rules for constructing /C„(t). 

COHERENCE FUNCTION FOR A CERTAIN CLASS OF INITIAL BATH STATES 

Eq. ((23)) determines the coherence function through the relation 

Ps.-+it) = Tp's'.Ut) = Vpi") (0)exp<' Ve^W } ■ (24) 



a a K n ^ 



We will now evaluate this for a certain class of initial bath states which is relevant for the electron spin decoherence 
problem. These states have the form p_b(0) = X^yPxx \x) {x\ (which corresponds to choosing n„ = \x) {x\)i where 
Ix) = ^f, \Ik,'ni^) is a product of cigenstates of the I^ (Ik{Ik + 1) and m^ are eigenvalues of I^ and /|). In this case, 
A/'q = 1. In other words, we arc focusing on initial bath states which are mixtures of states possessing well defined 
Ovcrhauser fields. 

The first step in evaluating Eq. ([24|l is to evaluate the correlators in Eq. ([23]) . 

TrB{{h+h-)'U^ih-h+y^-'^} = (xl {h-h+r^~\h+h-)'^ \x) . (25) 

For fixed values of qi and fc, this correlator is the sum of all possible ways of contracting raising operators I^ with 
lowering operators I^ . In general, a single contraction can involve an arbitrary number of operators (each sharing the 
same nuclear site index £), but in the large N limit, one can ignore all contractions involving more than two operators 
since these are far fewer in number relative to pairwisc contractions and should therefore have a subleading effect. 
This is tantamount to keeping only those contributions in which each raising and lowering operator pair has a site 
index which is distinct from every other pair. For example, one fourth-order contribution comes from the correlator 

{x\h-h+h-h+\x) = Y.il[AA {x\IijlliX\x)- (26) 

{^.} V J" / 

This correlator has two distinct pairwise contractions: ii = (-2, ^3 — ^4 and £1 = £4, £2 — £3- The first contraction 
gives the contribution 

Y: AlAl (xl lijl Ix) (xl liA IX) = E AlAl/-^^4f\ (27) 

where we have defined c\ = (x| Ii^^e |x), while the second contraction gives 

E AlAl (xl liJt Ix) (xl Ilji, Ix) = E ^'A'A^^^'^^'- (28) 

Therefore, the net result in the large N limit is 

(xl h-h^h-h^ Ix) = E ^\A\\^^t^T + 4f ^4^^"]- (29) 

To simplify this slightly further, we will invoke the large N limit once again to eliminate the inequality condition 
between £1 and £2- Then defining 

d± EE E Ahf^^ = E A'AWe + 1) - m}im} ± 1)], (30) 

€ e 

we may write 

{X\h-h+h-h+\x)^d+{d+ + d-). (31) 



A correlator involving qi pairs of raising and lowering operators will have qi ! different pairwise contractions which will 
each contribute to the coherence function. Each contraction will contribute a monomial involving dt and d~ , so that 
the full correlator will contribute some polynomial in these quantities. For example, the correlator with g^ = 3 and 
k = 2 involves six contractions and evaluates to 



(xl h-h+h+h-h+h- \x) = 2(4)2rf- + Ad+{d-f 



(32) 



When all of the correlators contributing to the coherence at a given order n are added together, the result is again a 
polynomial in d^ . Interestingly, this polynomial (for n > 2) is related to Eulerian numbers: 



E 



(-1) 



r+l r qi 



{q.}eV(n/2) 



nLi9.! rNj;, 



HE 



yTr.{(/.+.-)^-n„(/.-/.-r-}-(;^ 



fc=0 ^ ' 



(33) 

The symbol / ) denotes an Eulerian number, which for positive integer p and non-negative integer q can be expressed 
as 



9+1 



EM)^n7'^-^- + i 



i=0 



The sum over k in (|33p can be expressed as a polylogarithm: 
fc=0 ^ ' 



'^^-(4-rf^r/^E(|)^--^ 



(34) 



(35) 



The final expression involving the infinite sum is strictly speaking only valid for d7^ < d^ (which is true for negative 
polarizations '^lA'jrn^ < 0), but we will tacitly perform an analytic continuation at the end of the calculation. We 
may then express the xth coherence degree of freedom as 






exp 



4§:«+''- 



The sum over n is now easily done: 

fitjid+-d-)Y/^ 1 



j=i \ax/ ■? „>2 V 



2f7v 



fitj{d+ -d ) 



1 



V 2VLy / (n/2)! 

)i>2 ^ X / \ I ) 



exp 



itjid-^ ~d^ 



2n^ 



- 1- 



(n/2)! 

itj{d+ - d^) 
201 ' 



and the resulting sums over j are readily identified as logarithmic, so that 
^ ' I ~ 2^ \ or) ] ~ 



.d' 



E 

j=i 



1 ^ nt3{d+-d-) y/^ 

j^^[ 2n^ ) (n/2)! 



log 



4-^x 



dt -dx exp[2f|-(dj - d^)] 



2r2x '^ 



Adding this result to the second order (n = 2) contribution and exponentiating, we finally obtain 



pi"i+(o) 



d+-d^ 



'+„ 4n-^ 



d^e 



j-g4n^v-x- 



(36) 



(37) 



(38) 



(39) 



The full coherence function is obtained by summing over degenerate Overhauser states. To express the coherence 
function in the Schrodinger picture, we simply include the phase factor e*' + x'* inside the sum over x, where 
/if, = Tr^jIIp^,/!^}; this follows immediately from the property [11^, TJq] = 0. Denoting the Schrodingcr-picture density 
matrix as p, the Schrodinger-picture coherence function is then 



W{t) 



/5s,-+(0) 



E 



X "x s 



d^e'"- 



(40) 
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